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Universal molecules trapped with discrete scaling symmetries 



(N 

O 
(N 



00 



c/2 



O 

O 

(N 



X: 



Yusuke Nishida^ and Dean Lee^ 

^Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 
^Department of Physics, North Carolina State University, Raleigh, North Carolina 27695, USA 

(Dated: February 2012) 

When the scattering length is proportional to distance from the center of the system, two particles 
are shown to be trapped about the center. Furthermore, their spectrum exhibits discrete scale 
invariance whose scale factor is controlled by the slope of the scattering length. We also elucidate 
how the emergent discrete scaling symmetry is violated for more than two bosons, which may shed 
new light on Efimov physics. Our system thus serves as a tunable model system to investigate 
universal physics involving scale invariance, quantum anomaly, and renormalization group limit 
cycle, which are important in a broad range of quantum physics. 
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Introduction 

When particles attract by a short-range interaction 
with large scattering length, their low-energy physics be- 
comes universal [if. Ultracold atoms are ideal to study 
such universal physics because of the tunability of inter- 
atomic interactions and can provide insights applicable 
in a broad range of physics. One of the most striking 
phenomena in universal systems is the Efimov effect, i.e., 
formation of an infinite tower of three-body bound states 
characterized by discrete scale invariance (2|. Although 
the Efimov effect was originally predicted in the context 
of nuclear physics, it is now subject to extensive research 
in ultracold atoms Q. 

In this Letter, we propose novel systems in which two 
particles exhibit discrete scale invariance in their spec- 
trum. In order for this to happen, their interaction needs 
to be scale invariant Q . It is usually considered that the 
short-range interaction can be scale invariant only when 
the scattering length a is set to be zero or infinite. How- 
ever, there is another possibility: The scattering length 
is made space-dependent and tuned to be proportional 
to distance from the center of the system; a{x) = c\x\. 
This interaction is scale invariant because there is no 
dimensionful parameter in it. 

The emergence of the discrete scale invariance can be 
understood intuitively by using the Born-Oppenheimer 
approximation. Suppose one particle is much heavier 
than the other particle. With the heavy particle fixed at 
X, the light particle forms a bound state with binding en- 
ergy —h^/[2fj,a{x)'^], which in turn acts as an effective po- 
tential for the heavy particle. Therefore, one can design 
any attractive potential by tuning the space-dependence 
of the scattering length. In particular, when a{x) = c\x\, 
the effective potential becomes an inverse square poten- 
tial, for which it is well known that the spectrum exhibits 
discrete scale invariance. Since two particles are trapped 
about the center of the system with the discrete scaling 
symmetry, we shall call our system as a scaling trap. 

This conclusion can be established for any mass ratio 
by solving the two-body problem exactly with the space- 



dependent scattering length. While our idea works in any 
spatial dimensions, we shall give extensive and detailed 
analyses in one dimension and then present key results 
in two and three dimensions. 

Two particles in one dimension 

Two interacting particles in one dimension are de- 
scribed by the Schrodinger equation (hereafter ?i = 1): 



2M 2fi , 



^{X,x)^Ei;{X,x). (1) 



Here M = mi + TO2 and fi = mim2/(rni + 7712) are 
total and reduced masses and X — (mixi -\- m2X2)/M 
and X — xi — X2 are center-of-mass and relative co- 
ordinates. For a zero-range interaction whose strength 
depends on the position X, the interaction potential is 
written as V{X,x) — — ^J^x) ^i^)^ where a{X) is the 
space-dependent scattering length. For a bound state 
solution with E = — k^/(2M), the Schrodinger equation 
is formally solved by 
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where iplP^p) is the wave function in momentum space 
and = JdX e^"^^-^ Ijix) Fourier transform of 

the inverse scattering length. By integrating both sides 
of Eq. ([2|) over p, we obtain an integral equation solved 
byx(P)^/|f^(P,p): 
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We note that when a{X) = a > is uniform, there is a 
single bound state with binding energy \E\ — l/(2/ia^). 

Now for the linearly space-dependent scattering length 
l/a{X) = l/{c\X\), l/a(P) is ill defined because of 
the divergence at X = 0. This divergence needs to be 
regularized, for example, by a sharp cutoff l/a{X) = 
0{\X\ - e)/{c\X\) or by a smooth cutoff l/a{X) = 
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FIG. 1. Density distribution n±{xi) in the two-particle 
trapped state with even (left) or odd (right panel) parity in 
units of K — 1. Solid, dashed, dotted curves correspond to 
c = 0.5, 1, 1.5, respectively, with equal masses mi — m-2. 



l/(cVA^ + e^). In either case, the limit of an infinitesi- 
mal cutoff e leads to l/a(P) -(2/c) ln(e|P|), for 
which the analytic solution to Eq. ([3]) is obtained as 



sin[s±arcsinh(P/K) + 7r(l ± l)/4] 
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The upper (lower) sign corresponds to the even (odd) 

^±1 ^ 



parity channel and |iV±p = — — is the nor- 
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malization constant and s± solves 
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Note that ,s+ has a solution for any c > 0, while s_ has 
a solution only for < c < -I The latter range of c 

is assumed below unless otherwise stated. 

Then the inverse Fourier transform of x±(^) leads to 
the wave function with two particles at the same point: 
X±{^) = V'±(-^jO) cx Kis^{K\X\). An important obser- 
vation is that this wave function toward the origin X — 
oscillates as 7/'±(X, 0) — ?► cos[s± ln(K|X|/2) — argr(?s±)]. 
The phase of this oscillation is fixed by the precise be- 
havior of a(X) near X = 0, which is not universal and 
depends on experimental setups. However, what is uni- 
versal is that because of the logarithmic periodicity in k. 



if K = K-t is a solution, then k = e" 



K± are all solu- 



tions. Therefore, in each parity channel, there exists an 
infinite tower of two-body bound states characterized by 



discrete scale invariance: E. 
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This is exactly the same physics as the Efimov effect, 
while the difference should be emphasized that our bound 
state consisting of two particles is trapped about the cen- 
ter of the system (see Fig. [T] below) and the scale factor 
\± = e^l^^ is tunable by the slope of the scattering 
length as seen in Eq. ([5]). We also note that the full scale 
invariance demonstrated by the classical Hamiltonian ([T]) 
with a(X) = c\X\ is broken down to the discrete subset 
by the scale k± generated in quantum mechanics. This 
is known as a quantum anomaly. 

With the use of the wave function obtained from 
Eqs. the density distribution of a particle with 



FIG. 2. Momentum distribution PiP±(pi) in the two-particle 
trapped state. See Fig. [1] for details. 

mass TOi and its momentum distribution p±{pi) are plot- 
ted in Figs. [T] and [21 respectively, for c = 0.5, 1, 1.5 
with equal masses mi = m2. In particular, the momen- 
tum distribution has an oscillatory large momentum tail 
P±{Pi) K'^t±{pi) at \pi\/k oo given by 
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This logarithmic oscillation signals the discrete scale in- 
variance and exactly the same tail emerges in any few- 
body and many-body states as we will show later. In 
contrast to the Efimov effect in which the oscillatory tail 
appears at the subleading order @, Q , it appears at the 
leading order in our scaling trap. This will make its ob- 
servation easier by a timc-of-flight measurement in ultra- 
cold atom experiments. 

More than two particles 

A longstanding problem in Efimov physics is whether 
the discrete scale invariance demonstrated for three par- 
ticles persists for larger numbers of particles The 
scaling trap realizes a novel pattern of discrete scaling 
symmetry violation for bosons: Different particle sectors 
obey different scaling laws and incommensurate scalings 
among them result in the breakdown of discrete scale 
invariance. 

This pattern can be explained easily in the limit c 1 
where the Born-Oppenheimer approximation is applica- 
ble. Recall that one-dimensional bosons form an TV-body 
bound state with binding energy —N{N^ — l)/(()ma'^) 0. 
For a — c\X\, this binding energy acts as an effec- 
tive potential for the center-of-mass motion of the TV- 
body cluster. Accordingly, TV bosons form an infinite 
tower of trapped states characterized by a discrete scal- 
ing symmetry set by a scale factor \n — e^/"™ with 
SN = ^iV2(iV2 _ l)/(3c2) -t- 0(c°). However, these N- 
body trapped states for iV > 3 are actually unstable res- 
onances coupled with a continuous spectrum. For exam- 
ple, when = 3, there are continuum states composed 
of a free particle and two-body trapped state obeying a 
discrete scaling symmetry set by A2 . Because dilatations 
with respect to A2 ~ e'^'^^^ and A3 sa f^^'^/i'^^) are incom- 
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FIG. 3. A'^-boson resonance energy -E^' in units of jj. — 1 and 
unit lattice spacing versus excitation number n at c = 0.25. 
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mensurate, the discrete scaling symmetry breaks down 
for three bosons and hence more. 

Further insights into the nature of iV-boson resonances 
can be obtained in the same limit c <C 1. Because the 
local binding energy scales as — l/(c|X|)^, their relative 
wave function is localized within a separation ^ c\X\. On 
the other hand, their center-of-mass wave function is ^ 
Kispi{KN\X\), which oscillates rapidly at KAr|X| <g; 1 and 
decays exponentially at KAr|X| S> 1. This A^-body cluster 
can decay into a deeper TV'-body cluster with 2 < N' < 
N whose wave function also oscillates rapidly but with a 
different logarithmic period tt/sn' > tt/sn- Because of 
the resulting small overlap between their wave functions, 
A^-boson resonances are expected to have small decay 
widths and thus obey an approximate discrete scaling 
law set by Aat. 

To confirm this consideration, we numerically com- 
puted A^-boson resonance energies for A^ = 2, 3, 4. Here 
a Hamiltonian lattice formalism and iterative eigenvector 
methods were used in a semi-infinite system with a hard 
wall boundary at the origin. At c = 0.25, decay widths 
of resonances are indeed found to be negligible and scal- 
ing behaviors of resonance energies are clearly seen in 
Fig. [31 Their scale factors are extracted as A2 ~ 1.482(2), 
A3 « 1.172(2), and A4 ~ 1.10(2), which are in good 
agreement with the Born-Oppenheimer approximation; 

^ g7rc/V^'(^'-l)/3. 

In contrast to bosons, we did not observe any reso- 
nances in the spectrum of three fermions (two of one com- 
ponent and one of the other) with equal masses. This is 
indeed expected because two-component fermions with 
equal masses do not form any bound states with more 
than two fermions [7|. Therefore, their discrete scaling 
symmetry cannot be violated by the pattern elucidated 
above. 

On the other hand, there exists another pattern of 
discrete scaling symmetry violation which is common to 
bosons and fermions. For more than two particles, both 
of even and odd parity two-particle states contribute in 
general. Because they obey incommensurate scalings set 
by A4. and A_, the discrete scaling symmetry breaks 
down. However, this pattern does not take place in the 



range c > tt where only the even parity channel exhibits 
the discrete scale invariance [see Eq. ([5])] . Here it is pos- 
sible that an arbitrary number of fermions maintains the 
discrete scaling symmetry set by A+. Accordingly, we 
expect two-component fermions with equal masses to ex- 
hibit three phases with distinct symmetries as a function 
of the inverse slope; phases with full scale invariance at 



< 0, discrete scale invariance at < 



and no scale invariance at l/vr < . Their many- 
body physics and phase transitions in between will be 
extremely interesting and to be explored in the future. 

Effective field theory 

The scaling trap can be formulated in the language of 
effective field theories. A local field theory to be consid- 
ered is 

1 



H 



dx 



9+ 



9- 



^c\x 



(7a) 



-44^P2Mo) + ^v[^\4]^[^2^im, (7b) 

where i — 1,2 is summed and the space argument (x) 
acts on all operators on its left. The first part (iTal) de- 
scribes particles interacting by the zero-range interaction 
with the linearly space-dependent scattering length. As 
we discussed above, such an interaction is singular at 
the origin and needs to be regularized by introducing 
a cutoff. The independence of physical quantities from 
the cutoff is ensured by counter-terms. According to the 
above argument of discrete scaling symmetry violation, 
iV-body counter-terms, V'^^V'^(O) and V[V'^^]V[?A^](0), 
are needed for bosons with each N > 2. On the other 
hand, two counter-terms (I7bp are expected to be suffi- 
cient for an arbitrary number of two-component fermions 
with equal masses. 

An expression for the cutoff-dependent coupling g± (A) 
is obtained in the same spirit of Ref. 8]: We employ a 
sharp momentum cutoff \P\ < A and require that two- 
body physics becomes independent from the choice of A. 
The two-body sector of our field theory ([7]) is equivalent 
to the Schrodinger equation in Eq. ([T|) but with the in- 
teraction potential; V{X,x) = [-l/(c|A|) -I- g+S{X) + 
g-W x5{X)V x]5{x) / ^. Accordingly, the integral equa- 
tion in Eq. ([3]) is modified into that in which \/a{P—P') 
is replaced by -{2/ c)\n{\P - P'\/ K)- g+- g^PP' . Then 
we require that its solution x±{P) for |P| ^ A in Eq. @ 
does not change when the cutoff is changed from A to A'. 
This requirement is satisfied by choosing 

, . , 2 cot[0+(A)] , 

g+ (A) = a+ Lr+^ ^ 

c s+ 

.2 sin[0_(A)] + s_ cos[0_(A)] 
A g_(A) = q;_^— - — -— — --— , (8b) 

^ ' sin[(/)_(A)] - s_ cos[0-(A)] ^ ' 

where 0±(A) = s± \n{2K/ k±) + I3±. With numerical con- 
stants a± and /3± , the analytic expression (|5]) fits to a nu- 
merical solution accurately. An important observation is 
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that these couphngs run in logarithmically periodic ways 
as functions of the cutoff. This is known as a renormal- 
ization group limit cycle and the Efimov effect is its rare 
manifestation in physics. Our scaling trap is newly added 
to a short list of systems demonstrating the limit cycle. 

The field-theoretical formulation (O is useful to de- 
rive universal relationships valid in an arbitrary few- 
body and many-body state 0. An operator prod- 
uct expansion of J dxic-'P^'^^^liXi - ^)ipi{Xi + ^) 
at \pi\ — >■ OD is dominated by two lowest local op- 
erators; 0+{Xi) = VlV'lV'2V'i(^i) and O-(Xi) = 
'V[il;lipl]V[ip2'4'i]iXi). By matching their matrix ele- 
ments with respect to two-particle trapped states with 
even and odd parities, a Wilson coefficient of 0± is 
found to be -^^^t±{pi)6{Xi). Accordingly, the 
momentum distribution of a particle with mass mi ex- 
hibits the oscillatory large momentum tail; p{pi) — 
t-^.{pl)C+ + i_(pi)C_ in any state of the scaling trap. 
The form of each term is fixed by t± {pi ) obtained in 
Eq. ([6]) for equal masses, while its magnitude is set by 
a local contact density; C± = -^{-^^O±{0)). By 
applying the Hellmann-Feynman theorem to the Hamil- 
tonian ((T]), we find that C± measures how an energy of 
the state under consideration changes with respect to k,± : 
K±dE/dK± = -C±/M 9]. 

Scaling traps in two and three dimensions 

So far we have focused on physics in one dimension 
but our idea works equally in two and three dimensions. 
A two-body bound state problem in an arbitrary spatial 
dimension d reduces to solving an integral equation which 
is an analog of Eq. ([3]) in real space: 
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V'(X,0)/ln[Aa(X)] and x^{X) 



limi^i^n -!^\\x\ip{X ,x)\. The scattering length a{X) is 



Here X2{X) 

l|xKO 

defined so that when a{X) = a > is uniform, there is a 
single bound state with binding energy \E\ = l/(2/ia^). 

The emergence of a scaling trap for a{X) = c\X\ 
is deduced by whether there exist solutions that os- 
cillate as Xd{X) \X\^-'^^^' toward the origin 
K,\X\ — > 0. By substituting an ansatz X2{X) ^ 
\X\-^+''^^e'^^^ or X3(^) ^ |X|-2+*''«P^(cos6'^) into 

Eq. © with K — > 0, we find that si solves ln( ^) = 



Re[r'(^±i±H.)/r( , 



d ^ 2 



|P^£+2±Hi)/r(£+i±^)| in d = 3. For each angular 
momentum ^ where si has a solution, there exists an 
infinite tower of two-body trapped states characterized 
by discrete scale invariance: i?^"-* = — e~^"'^/*'«;^/(2Af). 
We expect similar patterns of discrete scaling symmetry 
violation for more than two particles as elucidated in one 
dimension. 



Remarks on experimental realization 

In this Letter, we proposed scaling traps in which two 
particles form an infinite tower of trapped states charac- 
terized by discrete scale invariance. The key idea is to 
make the scattering length proportional to distance from 
the center of the system. Such space-dependent interac- 
tions can be realized in ultracold atom experiments by 
spatially varying a magnetic- or optical-field intensity or 
by varying transverse confinement lengths along a longi- 
tudinal direction. 

If two particles correspond to different spin states of 
an fermionic atom, our two-body bound states with Efi- 
mov character are long-lived because three-body recom- 
binations are strongly suppressed by the Pauli exclusion 
principle Furthermore, a scale factor can be easily 
controlled by the slope of the scattering length. This will 
greatly facilitate an observation of the discrete scale in- 
variance, for example, by a radio-frequency spectroscopy 
or a time-of- flight measurement. Therefore, the scaling 
trap overcomes common difficulties in ultracold atom ex- 
periments of Efimov physics arising from an instability 
of three-body bound states and their sizable scale factor 
w 22.7. 

We also elucidated that the discrete scaling symme- 
try emergent for two particles is inevitably violated for 
three or more bosons by the appearance of resonance 
states while not for fermions. It is possible that insights 
developed here shed new light on Efimov physics. Our 
scaling trap thus serves as a tunable model system to 
investigate universal physics involving scale invariance, 
quantum anomaly, and renormalization group limit cy- 
cle, which are important in a broad range of quantum 
physics. 
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